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By  using  more  realistic  a  priori  knowledge  about  the  target  motion,  tracking  of 
maneuvering  targets  for  homing  missiles  is  enhanced.  Since  certain  targets  are 
assumed  to  execute  evasive  maneuvers  orthogonal  to  their  velocity  vector,  a  new 
stochastic  dynamic  target  model  is  proposed  where  this  orthogonality  is  embedded 
Along  with  this  new  acceleration  dynamic  model,  the  orthogonality  is  also 
enforced  by  the  addition  of  a  fictitious  auxiliary  measurement.  The  target 
states  are  estimated  by  the  modified  gain  extended  Kalman  filter  (MGEKF),  and 
the  angular  target  maneuver  rate  is  constructed  on-line.  A  guidance  law  that 
minimizes  a  quadratic  performance  index  subject  to  the  assumed  stochastic 
engagement  dynamics  that  includes  state  dependent  noise  is  derived.  This 
guidance  law  is  determined  in  closed  form  where  the  gains  are  an  explicit 
function  of  the  estimated  target  maneuver  rate  as  well  as  time  to  go.  The 
numerical  simulation  for  the  two-dimensional  angle-only  measurement  case 
indicates  that  the  proposed  target  model  with  the  MGEKF  leads  to  remarkable 
estimation  of  the  target  states.  Furthermore,  the  effect  on  terminal  miss 
distance  using  this  new  guidance  scheme  is  given. 
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Abstract 

By  using  more  realistic  a  priori  knowledge  about  the  target  motion,  tracking  of  ma¬ 
neuvering  targets  for  homing  missiles  is  enhanced.  Since  certain  targets  are  assumed 
to  execute  evasive  maneuvers  orthogonal  to  their  velocity  vector,  a  new  stochastic 
dynamic  target  model  is  proposed  where  this  orthogonality  is  embedded.  Along  with 
this  new  acceleration  dynamic  model,  the  orthogonality  is  also  enforced  by  the  ad¬ 
dition  of  a  fictitious  auxiliary  measurement.  The  target  states  are  estimated  by  the 
modified  gain  extended  Kalman  filter(MGEKF),  and  the  angular  target  maneuver 
rate  is  constructed  on-line.  A  guidance  law  that  minimizes  a  quadratic  performance 
index  subject  to  the  assumed  stochastic  engagement  dynamics  that  includes  state 
dependent  noise  is  derived.  This  guidance  law  is  determined  in  closed  form  where  the 
gains  are  an  explicit  function  of  the  estimated  target  maneuver  rate  as  well  as  time 
to  go.  The  numerical  simulation  for  the  two-dimensional  angle-only  measurement 
case  indicates  that  the  proposed  target  model  with  the  MGEKF  leads  to  remarkable 
estimation  of  the  target  states.  Furthermore,  the  effect  on  terminal  miss  distance 
using  this  new  guidance  scheme  is  given. 
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1.  Introduction 

The  target  tracking  problem  for  homing  missile  guidance  involves  the  estimation 
of  problems  of  estimating  large  and  rapidly  changing  target  accelerations.  The  time 
history  of  target  motion  is  inherently  a  jump  process  where  the  acceleration  levels 
and  switching  times  are  unknown  a  priori.  Due  to  this  arbitrary  and  unpredictable 
nature  of  target  maneuverability,  target  acceleration  cannot  easily  be  modeled. 

A  considerable  number  of  tracking  methods  for  maneuvering  targets  have  been 
proposed  and  developed  for  both  new  target  models  and  filtering  techniques[l]-[7]. 
In  spite  of  the  numerous  modeling  and  filtering  techniques  available,  target  accelera¬ 
tion  estimation  using  angle-only  measurements  is  relatively  poor.  Usually,  the  target 
tracking  problem  is  approached  by  modeling  target  acceleration  with  a  first-order 
Gauss-Markov  model  and  applying  the  extended  Kalman  filter(EKF).  One  difficulty 
with  the  Gauss-Markov  model  is  that  the  assumed  large  process  noise  spectral  den¬ 
sity  induces  Kalman  filter  divergence  even  when  the  target  maneuver  is  not  present 
and  can  lead  to  a  target  acceleration  magnitude  estimate  which  exceeds  the  actual 
maximum.  Another  problem  is  that  the  target  motion  is  not  well  represented  by 
Gauss-Markov  diffusion  process. 

In  an  effort  to  alleviate  these  problems,  the  circular  target  model  has  been  pro¬ 
posed  as  a  target  motion  model,  where  the  phase  angle  is  a  Brownian  motion  pro¬ 
cess  and  the  acceleration  magnitude  can  be  either  a  random  variable  or  a  bounded 
stochastic  process.  This  target  model  was  suggested  in  [4]  h-g  concepts  extracted 
from  [9,10]. 

By  including  a  priori  knowledge  of  the  target  motion,  improved  estimates  of  the 
target  states  can  be  obtained.  This  idea  is  included  in  [6,7]  by  using  a  target  acceler¬ 
ation  model  which  employs  a  target  mean  jerk  term.  For  conventional  targets  such  as 
winged  aircraft,  the  longitudinal  acceleration  component  is  often  negligible  compared 
to  the  lateral  component  in  evasive  maneuvers.  This  notion  fits  the  circular  target 
model  where  the  angular  rate  term  is  estimated  to  account  for  the  actual  dynamics  of 
the  coordinated  turn.  This  model  is  presented  in  Section  2.  However,  an  approximate 


state  expansion  is  required  to  handle  the  unknown  angular  rate  in  the  target  model. 
This  approximate  dynamical  system  used  for  estimation  is  presented  in  Sections  3.1 
and  3.2.  Furthermore,  in  Section  3.3  the  orthogonality  between  target  velocity  and 
acceleration  can  be  viewed  as  a  kinematic  constraint  where  compliance  is  inforced  by 
including  this  constraint  as  a  pseudo- measurement [7, 8].  The  approximate  target  dy¬ 
namics  and  pseudomeasurement  are  included  in  the  modified  gain  extended  Kalman 
filter(MGEKF)[ll]  and  is  presented  in  Section  4.  The  MGEKF  is  selected  because 
of  its  superior  performance  over  the  EKF  especially  for  bearing-only  problems.  In 
Section  5,  a  linear  quadratic  guidance  law  is  derived  for  this  circular  target  model. 
This  guidance  law  remains  a  linear  function  of  the  estimated  states,  but  the  guidance 
gains  obtained  in  closed  form  are  a  nonlinear  function  of  the  estimated  rotation  rate 
and  time  to  go.  Finally,  a  numerical  simulation  is  performed  for  a  two-dimensional 
homing  missile  intercept  problem.  Both  the  estimation  process  and  the  terminal  miss 
are  enhanced  by  the  new  models  and  the  associated  estimator  and  guidance  law. 

2.  Target  Acceleration  Model 

In  this  section  the  circular  target  acceleration  model  is  presented,  and  the  dynamic 
consistency  between  this  target  model  and  an  assumed  nonlinear  target  model  is 
discussed. 

2.1.  Circular  Target  Motion 

The  two-dimensional  homing  missile  guidance  scenario  is  described  by  two  sets  of 
nonlinear  dynamic  equations  of  motion  for  the  missile  and  target 

xm  =  Vmcos9m  ,  ij  = 

ilM  =  VMsin6M  ,  yT  = 

Vm  =  am,  ,  Vt  = 

,  0T  = 

where  (zy  ,  ip)  and  (yM  ,  Vt)  are  inertial  coordinates,  V\j  and  Vj  are  the  velocities, 
aw,  and  ar  are  the  accelerations,  and  9m  and  9j  are  the  flight  path  angles 
[Fig.  2].  The  subscripts  M  and  T  denote  the  missile  and  target,  respectively,  and 
aM,  and  dMn  are  tangential  and  normal  accelerations,  respectively.  Only  the  normal 


Vjcos9t 

Vjsin9j 

0 

dj'/Vj' 


(1) 
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component  of  the  acceleration  contributes  to  changing  angular  orientations  of  each 
vehicle,  and  the  target  is  assumed  to  fly  at  constant  speed. 

The  following  target  model  is  assumed  to  be  used  in  the  filter.  The  objective  is 
to  choose  a  model  that  is  linear  in  order  to  reduce  the  numerical  computation  of  the 
filter,  but  reasonably  consistent  with  the  nonlinear  model  so  that  the  estimates  are 
of  good  quality.  The  target  model  for  the  filter  in  two  dimensions  is 

aTs(t)  =  aTCos(ut  +  0),  dTt(t)  =  arsin^ut  +  0)  (2) 


where  ay  is  a  constant  which  is  unknown  a  priori,  ui  is  the  angular  velocity  to  be 
estimated  in  a  right-handed  coordinate  system,  and  0  is  a  Brownian  motion  process 
with  statistics 

E[d0\=  0,  E[dO*]  =  Qdt,  0  =  1/r©  (3) 

Here,  0  is  the  power  spectral  density  of  the  process  and  t©  is  the  coherence  time, 
the  time  for  the  standard  deviation  of  0  to  reach  one  radian.  While  in  the  previous 
circular  target  model[4]  the  acceleration  components  were  just  a  diffusion  process 
along  a  circle,  those  in  the  new  model  are  related  to  the  actual  target  motion  through 
a  term  of  physical  meaning,  u. 

2.2.  Dynamic  Consistency 

In  order  to  see  how  the  current  model  approximates  the  assumed  nonlinear  target 
dynamics(l),  consider  a  deterministically  equivalent  case.  Integration  of  (2)  with 
0  =  0  and  u>  >  0  yields 

Vx  =  —sinut,  Vv  =  ——cosuit  —  Vt  +  — .  (4) 

u  ui  u 

where  the  initial  conditions  are  Vs  =  0  and  Vy  =  —  Vp.  Adding  the  square  of  each 
component  and  moving  all  terms  to  the  left  hand  side  gives 

2(1  -  cosut)l(— )2  -  —  VT]  =  0.  (5) 

U)  Ul 

For  this  equation  to  hold  for  all  t  >  0 


u  = 


aj 

Vt 


(6) 
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which  is  equivalent  to  the  differential  equation  for  angular  rate  in  (1). 

Furthermore,  by  taking  the  dot  product  of  the  velocity  vector  with  the  acceleration 
vector,  we  obtain 

Vt-St  =  otI—Vt  +  — 1  sinut  =  0  (7) 

using  (6)  for  all  t  >  0.  This  orthogonality  between  target  velocity  and  acceleration 
demonstrates  the  dynamic  consistency  of  the  proposed  target  acceleration  model  for 
the  filter  with  the  nonlinear  target  dynamics. 

3.  New  Dynamic  and  Measurement  Models  for  Estimation 

The  previous  section  dealt  with  a  new  circular  model  for  filter  implementation  in 
order  to  exploit  an  assumed  characterization  of  the  motion  of  a  typical  target.  In 
this  section,  the  stochastic  dynamic  equations  for  the  new  target  model  are  derived. 
Furthermore,  the  kinematic  fictitious  measurement  suggested  in  [7,8]  is  also  discussed. 
3.1.  Formulation  and  Approximation 

Ito  stochastic  calculus[12]  applied  to  the  Eq.(2)  results  in  a  stochastic  differential 
equation  with  white  state-dependent  noise[9]. 


dar z 

•-?  ->• 

ar. 

'  0 

-dO  ' 

aTl 

= 

dt  -f- 

darv 

.  -  -f. 

.  aTy  . 

dQ 

0 

.  aTy 

where  the  elements  —  f  in  the  drift  coefficient  are  the  Ito  correction  terms.  Note 
that  the  problem  is  nonlinear  due  to  the  unknown  u>.  To  avoid  solving  the  nonlinear 
problem,  u>  is  approximately  removed  by  an  expansion  of  state  variables  as  in  [10]. 
Define  new  states  as 


uaT 

>  al 

=  , 

t OCLy 

,  a] 

=  ua\  , 

(9) 


with  the  assumption 

M  =  « i.  (io) 

a 1 

By  augmenting  the  dynamics  of  these  new  states  to  (2),  an  approximate  dynamical 
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model  which  includes  this  new  target  model  is  truncated  as 


Note  that  u>  does  not  appear  explicitly  in  (11),  and  (11)  is  a  linear  stochastic  dif¬ 
ferential  equation.  An  idea  of  this  sort,  given  in  [10]  for  a  scalar  problem,  led  to 
significantly  improved  filter  performance. 

3.2.  Two-dimensional  Intercept 

The  two  dimensional  intercept  problem  is  developed  in  a  relative  inertial  coordi¬ 
nate  system.  The  system  dynamics  are  expressed  in  the  following  set  of  equations: 


=  «r,  yr  =  ur, 

«r  =  aTir  -  aMir ,  tv  =  aTyr  -  aMyr ,  (12) 

aTXr  =  a.TCOs(ujt  -f  9),  cLTyr  =  arsin(ut  4-  9). 

With  the  assumption  of  u  1  and  truncation  of  the  target  dynamics  up  to  the 
second  order,  the  ten-element  state  vector  is  defined  as 


1  =  (  Sr  yr  «r  vr  ax  ay  a\  a}y  a2y  ]T 

=  [  Xi  X2  X3  X4  XS  X6  Xj  Ig  x9  x10  ]T. 


Thus,  in  terms  of  the  expanded  state  space,  the  linear,  stochastic  state  differential 
equation  is  described  by 

dx  =  (Fx  +  Bu)dt  +  Gd9  (14) 

where  F  is  given  by  the  relations  in  (12)  and  by  the  coefficient  matrix  of  x  in  (11) 
where  i  =  2,  B  is  a  10  x2  matrix  of  zero  except  for  B3X  =  B42  =  — 1,  u  =  {a^tr ,  aMVr)T, 


G  =  [  0  0  0  0  -  x6  x5  -  x8  x7  -  xjo  x9  ]T. 


Measurement 


6 


Angle  measurement 

Angle  information  in  discrete-time  is  assumed.  Then,  the  measurement  at  time  tk 
is 

zi(tk)  =  hi(x(tk))  +  vk  =  tan~\yr/xr)  +  vk  (16) 

where  vk  is  a  white  random  sequence  with  statistics 

E[vk]  =  0,  £[u*UfT]  =  Vx8kl.  (17) 


Fictitious  measurement 

In  Ref. [7 ,8]  the  filter  performance  is  improved  by  introducing  a  kinematic  con¬ 
straint  based  on  a  priori  knowledge,  which  is  implemented  in  the  form  of  an  aug¬ 
mented  fictitious  measurement.  In  particular,  the  acceleration  vector  is  assumed  to 
be  related  to  the  velocity  as 

Vt  •  a-p  =  0  (16) 

under  the  assumption  that  the  target  accelerates  predominantly  orthogonally  to  its 
velocity  vector.  When  this  condition  is  not  met,  the  acceleration  has  a  component  in 
the  direction  of  the  target  velocity  as 

Vt  •  a  t  =  t]  (19) 

where  Vt  and  dr  are  assumed  to  be  random  vectors  representing  target  velocity 
and  acceleration,  and  rj  is  the  uncertainty  in  the  orthogonality.  This  idea  can  be 
implemented  in  the  form  of  a  discrete  pseudo-measurement  as 

z2{tk)  =  h2{x(tk))  +  T]k 

=  Tr,r  aTXr  +  VrVraTVr  +r]k  (20) 

where  Tjk  is  a  white  random  sequence  with  assumed  statistics 

E[r]k)  =  0,  E[rjkr1J]  =  V28kl.  (21) 
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Note  that  the  variance  of  the  measurement  noise  corresponds  to  the  tightness  of 
the  constraint.  In  other  words,  the  larger  the  variance,  the  more  relaxed  are  the 
requirements  on  longitudinal  acceleration.  This  fictitious  measurement  is  used  along 
with  the  angle  measurement  in  the  modified  gain  extended  Kalman  filter  described 
in  the  next  section. 

4.  Estimation  of  Target  States 

In  this  section,  the  modified  gain  extended  Kalman  filter(MGEKF)[ll]  is  derived 
for  the  circular  target  model  and  for  the  fictitious  and  angle  measurement  defined  in 
the  previous  section.  Also,  considered  is  the  method  to  reconstruct  the  maneuver  rate 
using  the  estimated  states.  Given  the  continuous-time  dynamics  and  discrete-time 
measurements,  as  in  the  previous  section,  construction  of  the  filter  is  completed  by 
specifying  the  time  propagation  and  measurement  update  procedures. 

4.1  Time  Propagation 

The  state  estimate  ir(i/*l_1)is  propagated  from  the  current  time  j  to  the  next 
sample  time  t,  by  integrating 

Ht/ti- 1)  =  Fx(t/ti _0  +  Bu{t), 

Ht/ti-i)  =  FP(t/ti_i)  +  P(t/ti_x)FT  +  E[GQGT],  (22) 

X(t)  =  FX(t)  +  X(t)FT  +  E[GQ(^}. 

given  the  posteriori  estimate  =  x(f,_i)  and  posteriori  pseudo-error  vari¬ 
ance  jP(t,_i/t,_i)  =  The  notation  R(t/ti_x)  denotes  the  value  of  some  quan¬ 

tity  R  at  time  t  given  the  measurement  sequence  up  to  time  The  integration  of 
the  covariance  of  the  state  X(t)  begins  with  X(0)  at  time  t  =  0.  Upon  integrating 
the  equations  above  to  the  next  sample  time,  the  propagated  estimates  are  obtained 
as  follows: 

x(ti)  =  P(u)  =  P(ti/t ,_i)  (23) 

It  should  be  noted  that  since  the  process  noise  is  state- dependent  ,  the  integration  to 
propagate  P(t)  also  requires  the  integration  of  X(t),  where  E[GQGT]  matrix  turns 
out  to  have  nonzero  elements  for  its  lower-right  6x6  matrix.  Note  that  the  approxi¬ 
mation  technique  reduces  the  originally  nonlinear  dynamics  to  linear  dynamics.  This 
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allows  for  the  closed  form  solutions  of  the  propagation  of  the  estimates  rather  than 
performing  on-line  integration. 

4.2  Measurement  Update)!!] 

The  states  axe  updated  as  follows: 


x(u)  =  + 

KM  =  P(t,)H(tt)T{H(U)P(t,)HT(i,)  +  V]-' 


(24) 


where 


V  = 


Vx  0 

0  V2 


and 


(25) 


#11, 

0, 


with 
#n  = 


Vr 


#12  = 


#12,  0, 

0,  #23, 

X, 


H, 


0, 

24, 


,0, 

#25, 


#■ 


26, 


_X?  +  y2  _ 

#23  =  X5,  #24  =  X6,  #25  =  13+14,,, 


#26  =  X4  +  V, 


y  m 


where  the  missile  acceleration  in  the  x  and  y  directions,  a\fz  and  aMy,  are  assumed 
to  be  measured  very  accurately  with  on-board  sensors. 

The  measurement  update  of  the  pseudo-error  variance  is  performed  by 


P(U)  =  [I-  K(U)g(z(U ),  x(u  ))}P(U)[I  -  K(U)g(z(U ),  x(i,))]r 
+K{ti)R(ti)K{ti)T 


(26) 


where  <7(2(4),  x(<,))  is  used  in  the  update  of  P  rather  than  #  of  Eq.  (25)  and  is  given 
as 


h(x(t,))  -  h(x(t,)) 


t  _-xVt{U)  1  _1  J/r(4) 

tan  — —  —  tan  _  — 
xr(t, )  xr(t.) 

Mx(t,'))  -  h2K*{ti)) 


=  sr(2(4),x(4))(x(4) -x(4))  (27) 


Note  that  g  is  a  2  x  10  matrix  of  function  explicit  only  in  the  known  quantities  z 
and  x.  In  this  sense,  the  function  h  has  a  universal  linearization  with  respect  to 
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the  measurement  function  z.  Unfortunately,  this  type  of  linearization  with  respect 
to  the  measurements  occurs  for  only  a  few  functions.  It  is  applicable  to  angle  mea- 
surements[ll],  but  not  for  our  new  pseudo-measurement.  Therefore,  we  must  for  the 
pseudo-measurement  revert  back  to  the  extended  Kalman  filter  form  and  define 


#(*&),  *(*0)  =  A2(x(*i))  -  h2(r{ti)) 

-  ah'  W*<<0  -  *((<)) 


dhi 


where  the  expression  for  U=g  *s  found  in  the  second  row  of  H  in  (25) 
For  angle  measurements^  1] 

~  =  9i{z{U),  x(t,))(x(t,)  -  x(f,)) 

=  -E{U)H{U){x{U)-x{U)) 


(28) 


(29) 


where 


m  = 

D{U)  = 
a(t.)  = 
H(z(U))  = 


D{ti)tan~xg(ti ) 

<*(*«)  ~ 

'fxr{U)2  +  yT(t,)2 
(xr(f,)xr(t,-)  4-  yr{ti)yT(ti)) 

DiumtiXz-itiMti) 

[sin  z(i<),  ,  -cos  z(t,)  ,0,0,0, ...] 


(30) 


As  discussed  in  [11],  g(z(t,),x(tj))  is  only  used  in  the  update  of  P(ti)  but  not  in  the 
gain,  since  it  was  empirically  shown  that  this  procedure  leads  to  an  unbiased  estimate 
of  the  state. 


4.3.  Estimation  of  u>  and  T, 
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Since  the  target  angular  velocity  term  is  embedded  in  the  states,  u ;  should  be 
reconstructed  using  the  estimated  states.  A  simple  way  to  determine  the  value  of  u> 
is  to  divide  the  states  as  w  =  —  or  ~.  However,  since  the  expanded  state  space  is 
originally  an  approximated  state  space,  this  might  lead  to  numerical  errors,  especially 
when  the  higher  approximated  terms  are  used.  By  relying  on  the  definition  of  the 
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vector  relation  between  and  velocity  and  acceleration,  the  target  angular  velocity 
can  be  obtained  without  using  the  augmented  states  for  approximation.  From  the 
assumed  dynamics  the  target  angular  velocity  during  its  evasive  maneuver  is 


Vy  x  Sy 


Thus,  by  using  the  state  estimates,  ui  is  constructed  as 


(31) 


Cj  =  sign(Vyxayy  -  Vr,dr,)| 


ay  | 

w 


(32) 


To  be  used  later  in  the  controller,  an  estimate  of  time-to-go,  Tg0, 
approximated  here  as 


A  A 

JL  _  R 

\R\  ~  I*  •  v\/R 


is  required,  and 


(33) 


where  R  and  R  are  the  estimates  of  relative  range  and  range  rate,  respectively,  and 
the  vectors  X  and  V  are  the  estimates  of  relative  position  and  velocity,  respectively. 


5.  Linear  Quadratic  Guidance  Law 

Based  on  the  estimated  states  and  the  estimate  of  the  rotation  rate  constructed 
from  the  estimated  states,  a  guidance  law  can  be  mechanized.  In  the  following,  a 
stochastic  guidance  law  is  determined  which  minimizes  a  quadratic  performance  index 
subject  to  the  stochastic  engagement  dynamics  including  the  stochastic  circular  target 
model(8)  under  the  assumption  that  states  including  the  target  states  and  the  target 
rotation  rate  are  known  perfectly.  This  assumption  simplifies  the  derivation  of  the 
guidance  law  enormously,  and  for  this  homing  problem  it  is  shown  that  the  solution  to 
the  stochastic  control  problem  with  state-dependent  noise  can  be  obtained  in  closed 
form.  The  solutio  1  Stained  does  not  produce  a  certainty  equivalence  controller  since 
the  guidance  law  explicitely  depends  upon  the  system  statistics. 

Note  that  since  the  noise  in  each  cartesian  direction  is  correlated  in  the  stochastic 
circular  target  model(8),  and  that  with  ayx  and  ayy  dynamically  coupled  through  u> 
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term,  the  guidance  commands  in  the  x  and  y  direction  cannot  be  achieved  indepen¬ 
dently.  Thus,  the  optimal  stochastic  controller  for  circular  target  model  is  based  on 
the  minimization  of  the  performance  index 


J  =  + »?]  +  f  Jl'K.  + 

subject  to  the  following  stochastic  system  of  linear  dynamic  equations 


dx 

dy 

du 

dv 

dar,, 

dary 


udt 

vdt 

(<Vr,  —a\f,)dt 

(ar  -  <tMy)dt 

(-|aT,  -uaTy)dt  -  aTydd 

(—%aTy  +  uarM)dt  -f  aixd0 


(34) 


(35) 


where  9  is  a  Brownian  motion  defined  earlier  and  £[•]  stands  for  an  expectation 
operator.  In  the  construction  of  filter,  the  inherent  nonlinearity  of  the  target  model 
was  removed  by  an  expansion  of  state  variables.  However,  for  the  guidance  law 
formulation,  the  rotation  rate,  w,  is  assumed  known,  although  it  must  be  constructed 
on-line  from  the  state  estimator(32). 

For  brevity  of  notation,  define  the  state  and  control  vectors  as  follows  : 


x  =  [x,y,u,v,aTl,aTy]T 

«  =  K,aMyf  (36) 

Then,  the  stochastic  control  problem  is  to  find  u  which  minimizes 

J  =  E{\Jq  uTRudT  +  \xTjSsx}}  (37) 

subject  to  the  stochastic  differential  equation  with  state  dependent  noise 

dx  =  [Ax  +  Bu]dt  +  D(x)d9  (38) 
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where 
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where  c  >  0.  Note  that 

D(x)  =  j2xjDj,  Dj  6 

j=i 

where  Xj  is  the  jth  element  of  x  and  where 


Du  Z?3,  D\ 


'  0  ' 

■  0  ' 

•  0  ■ 

0 

0 

0 

0 

0 

,Z?5  = 

0 

0 

,A>  = 

0 

0 

0 

0 

-1 

0 

.  1 . 

0 

(39) 


(40) 


(41) 


To  obtain  an  optimal  control  for  this  class  of  problem,  dynamic  programming[13,14] 
is  employed  where  the  Hamilton-Jacobi-Bellman  equation  becomes 


0  =  J°(x ,  t)  +  Min{Jx(Ax  +  Bu )  +  \[xT  A(JXXi  t)x  +  uTRu ]}  (42) 

u  2 

where  J°  is  the  optimal  return  function  and  the  subscripts  denote  partial  derivatives. 
The  elements  of  the  matrix  A  for  any  symmetric  matrix  W  is  defined  as 


A  =  tr[Di(t)TWDj{t)} 


(43) 


The  minimization  operator  in  (25)  produces 


u  =  -R~lBTJ°x 


(44) 
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(45) 


By  substituting  (44)  into  (42),  the  dynamic  programming  equation  becomes 

0  =  J°t  +  J°xAx  -  ±J°XBR-'BJ°X  +  i*TA (J°xx>Q,t)x 

The  optimization  problem  is  solved  by  explicitely  showing  that  the  equation  above 
has  a  solution  .  Asume  J°(x,t)  =  | xTS(t)x ,  then 

Jt  =  ±xTSx,  J°x  =  xTS,  J°XX  =  S  (46) 

With  this  assumption,  the  dynamic  programming  equation  is  satisfied  for  all  ?€  Rn 
if 

S  +  SA  +  AtS+A-SBR~1BtS  =  0,  S(tj)  =  Sf  (47) 

The  desired  optimal  controller  becomes 

u  =  —R~1BtSx  (48) 


where  S  is  the  solution  of  the  Riccati  equation  and  the  A(S,t)ij  =  tr(DfSDj\  leads 
to 

■  zero 

A  =  elements  (49) 

<i>66  “*^S6 

—  $56  SsS  . 

The  fact  that  A  has  only  nonzero  elements  for  its  lower-right  2x2  matrix  allows 
a  tractable  closed-form  solution.  To  see  the  characteristics  of  the  solution  in  a  simple 
manner,  matrices  are  partition  such  that  their  lower-right  block  partitioned  is  a  2  x  2 
matrix.  Then 
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where  the  block  matrices  satisfies  the  decomposed  Kiccati  equation 
—Su  =  ■S'n-^n  +  —  SiiBiR^Bi  Sn 

—S12  =  S’u/lij  +  5’i2-<422  +  -^n*?i2  —  SnBiR"1  Bf  S\2  (52) 

—^22  =  ^12^12  +  S22A22  +  ^22*^22  +  ^12*^22  "b  S  —  SuB\R  1 Si 2 

Since  the  S22  block  does  not  affect  the  block  matrices  Sn  and  Su,  the  optimal  control 

law  is  not  dependent  on  £22  •  Therefore,  the  closed  form  optimal  guidance  law  for 
this  special  class  of  problem  can  be  obtained  by  integrating  the  Riccati  equation 
backwards  without  requiring  the  explicit  evaluation  of  the  A  term.  In  particular,  the 
stochastic  optimal  control  problem  essentially  degenerates  to  a  deterministic  optimal 
control  problem  although  the  Ito  terms  are  retained.  The  solution  process  for  this 
deterministic  control  problem,  explained  in  detail  in  the  Appendix  A,  produces  a 
guidance  law  in  closed  form.  Note  that  the  deterministic  coefficient  A22  includes  the 
statistic  0.  Therefore,  the  resulting  controller  is  not  a  certainty  equivalence  controller. 
The  new  controller  becomes 


x{t) 

y(t) 

“(0 

v(t) 

aTx  (0 
L  °ry(  0 

where  the  gains  cx  to  c4  are  an  explicit  function  of  Tgo,  u,  and  0  as 


Cl 

0 

c2 

0 

c3  C4 

aMy 

0 

Cl 

0 

C2 

C4  C3 
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c  =  c-  {22kefr„  _  ainuTi'_g___ 


u>Q  .  ,  _  e-r..,  - “>') 


c\  =  c*  {-c dTgoe7T'°  +  [e^T,°  —  coswTgo) 


+  (cosu,Tgo-e*T'°)^~^}  (54) 
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Fig.  1  is  a  block  diagram  for  an  adaptive  guidance  scheme  for  a  homing  missile. 
Note  that  guidance  gains  are  functions  of  fg0 ,  estimated  time  to  go,  the  statistic  0  and 
the  estimated  maneuver  rate  u>.  Therefore,  for  the  bearing-only  measurement  system 
although  the  resulting  stochastic  guidance  law  is  sub-optimal  since  the  measurements 
are  nonlinear  functions  of  the  states,  the  explicit  dependence  on  the  estimate  of  the 
target  maneuver  rate  is  a  new  feature  which  should  help  reduce  terminal  miss  distance. 

6.  Numerical  Simulation 

For  a  particular  engagement  scenario,  the  performance  of  the  estimator  using  the 
new  target  models  and  that  of  the  guidance  law  are  evaluated. 

6.1.  Missile  and  Target  model 

Both  target  and  missile  are  treated  as  point  masses  and  are  considered  in  two- 
dimensional  reference  frames  as  shown  in  Fig.  2.  The  missile  represents  a  highly 
maneuverable,  short  range  air-to-air  missile  with  a  maximum  normal  acceleration  of 
lOOg's.  It  is  launched  with  a  velocity  M  =  0.9  at  a  10,000 ft  altitude  with  zero 
normal  acceleration.  After  a  0.4  sec  delay  to  clear  the  launch  rail,  it  flies  by  the 
guidance  command  provided  by  the  linear  quadratic  guidance  law  of  Section  5.  Also, 
to  compensate  for  the  aerodynamic  drag  and  propulsion,  the  missile  is  modeled  to  have 
a  known  longitudinal  acceleration  profile  :  a\f  =  25g's  for  t  <  2.6 sec,  a\f  =  —15 g's 
for  t  >  2.6sec.  The  target  model  flies  at  a  constant  speed  of  M  =  0.9,  and  at  an 
altitude  of  10,000/t.  It  accelerates  at  9  g's  either  at  the  beginning  or  in  the  middle 
of  the  engagement. 

The  actual  states  are  first  used  in  the  guidance  law  to  produce  consistency  in 
evaluating  the  performance  of  the  filter.  In  evaluating  the  actual  miss  distances,  the 
filter  state  estimates  are  used  in  the  guidance  law.  Two  engagements,  considered  in 
the  following  section,  are  shown  in  Fig.  3.  With  R\  and  R\f  denoting  initial  range  and 
maneuver  onset  range,  respectively,  engagement  1  is  the  situation  where  the  target 
maneuver  starts  at  the  beginning,  and  for  engagement  2  the  maneuver  starts  in  the 
middle. 
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Integration  of  actual  trajectories  is  performed  by  a  fourth-order  Runge-Kutta 
integrator  with  step  size  0.02  seconds.  The  variance  for  the  angle  measurement  is 
chosen,  as  given  in  [4],  to  be 

Vx  =  a  Va,  V0  =  +  5.625  *  10-7)/  At  rad2  (55) 

where  R  is  range,  At  is  filter  sample  time,- and  a  is  parameter  which  is  used  in  the 
simulation  indicating  different  levels  of  sensor  accuracy. 

As  mentioned  earlier,  the  variance  for  the  pseudo- measurement  can  be  interpreted 
to  show  how  strictly  the  orthogonality  assumption  between  the  target  velocity  and 
acceleration  is  to  be  kept.  By  allowing  some  acceleration  in  the  longitudinal  direction, 
a  reasonable  estimate  of  the  variance  to  be  used  can  be  given.  Suppose  that  the 
acceleration  component  in  the  velocity  direction  has  a  normal  distribution  with  zero 
mean.  Then  with  probability  0.95,  a  1  g  acceleration  while  flying  with  VT  =  970 ft/ sec 
leads  to  2  er  =  3.12  *  104  [ft2 /sec?],  where  a  is  the  standard  deviation,  which  results 
in  a  variance  V2  =  2.44  *  10s  [ft2 /sec?]2. 

Unless  otherwise  stated  the  filter  is  initialized  at  launch  with  the  true  relative 
position  and  relative  velocity  component  values.  Hence,  the  initial  values  for  the 
diagonal  elements  of  the  covariance  matrix  associated  with  position  and  velocity, 
P\l,  P‘2‘2,  P 33i  and  P44  are  set  to  ten.  On  the  other  hand,  little  knowledge  about 
target  acceleration  is  assumed  to  be  provided  at  the  beginning.  Therefore,  the  initial 
values  for  the  target  acceleration  and  expanded  states  are  zero.  Initial  values  of  the 
covariance  matrix  associated  with  target  acceleration  is  calculated  by  resorting  to 
the  definition  of  the  target  acceleration  at  t  =  0  given  in  (2).  Those  covariances 
are  produced  in  the  Appendix  B.  The  target  is  expected  to  execute  a  maximum 
acceleration  turn  in  its  evasive  motion,  and  the  missile  has  no  knowledge  about  the 
direction  of  target  rotation.  Note  that  0  is  a  Brownian  motion  process  beginning  at 
5(0)  =  9,  the  expected  angle  the  target  acceleration  vector  makes  with  respect  to  the 
x,  axis  at  the  time  of  launch,  and  arm„  is  the  expected  maximum  acceleration  of 
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target.  For  the  simulation  with  a  coordinate  system  having  one  axis  perpendicular  to 
the  initial  Vt  direction,  9  is  zero.  Then,  the  possible  nonzero  elements  of  the  upper 
triangular  part  of  the  initial  covariance  matrix  are  Pss(O),  Ps7(0),  Ps9(0),  Prr(O), 
^79(0),  and  ^99(0).  Furthermore,  no  information  is  available  about  the  direction  of 
maneuver,  and  the  possible  maximum  rotation  rate  can  be  either  positive  or  negative. 
Thus,  the  odd  powers  of  u>  are  taken  as  zero.  This  leaves  only  Pss(O),  Ps9(0),  Pn(0), 
Pgg(O)  as  the  nonzero  elements.  However,  a-value  of  ten  is  assigned  to  Pea,  Pea,  and 
P10 10  to  ensure  positive  definiteness  of  the  covariance  matrix  at  the  initial  time. 

6.4.  Filter  Results 

The  results  in  this  section  are  the  product  of  a  Monte  Carlo  analysis  consisting  of 
ten  filter  runs.  Along  with  the  miss  distance  calculations,  the  plots  of  the  estimation 
error  and  the  u>  estimates  versus  time  are  mainly  considered.  The  errors  are  calculated 
as  [P[ee]2  +  E[tvY\xP  where  E[tx\  and  E[ey]  are  the  averaged  values  of  errors  over 
ten  simulation  runs. 

Figs.  4-5  represents  the  results  for  the  engagement  1  where  the  target  maneuver 
initiates  at  t  =  0  and  the  pseudo-measurement  is  not  used.  As  seen  in  Fig.  4,  when 
there  is  no  switching  of  the  direction  of  the  target  acceleration  during  the  maneuver, 
the  estimation  results  get  better  as  the  power  spectral  density  of  the  process  noise 
decreases.  When  0  is  relatively  large,  the  target  acceleration  estimation  is  poor,  and 
divergence  of  the  position  and  velocity  estimation  occurs.  This  causes  the  estimate 
of  u )  to  deteriorate  as  time  goes  on.  However,  with  small  0,  the  filter  performance 
improves  considerably.  As  shown  in  Fig.  5,  estimation  improves  with  better  angle 
measurements. 

When  the  auxiliary  pseudo- measurement  is  also  implemented  in  the  filter,  estima¬ 
tion  performance  improves  over  the  case  when  only  an  angle  measurement  is  used. 
This  is  shown  in  Fig.  6  where  again  the  target  starts  its  acceleration  maneuver  at 
the  beginning  of  the  engagements  =  Rm).  At  first,  the  filter  with  the  fictitious 
measurement  seems  to  work  a  little  worse  than  the  filter  with  angle-only  measure¬ 
ment.  Then,  the  fictitious  measurement  promptly  works  as  if  it  suppressed  or  delayed 


18 


the  filter  divergence.  Note  that  the  effect  of  two  values  of  pseudo-noise  variance  are 
shown. 

The  role  of  the  fictious  measurement  is  more  observable  for  engagement  2  where 
the  target  maneuver  begins  in  the  middle  of  the  angagement(i?^  =  4000 ft).  As 
plotted  in  Fig.  7,  the  filter  equipped  with  only  the  angle  measurement  diverges  as 
soon  as  the  target  maneuver  occurs.  On  the  other  hand,  when  the  filter  is  augmented 
with  the  fictitious  measurement,  it  works  very  effectively.  The  divergence  of  position 
and  velocity  is  noticeably  suppressed,  and  the  acceleration  estimate  tend  to  return  to 
its  actual  value  from  an  instantaneous  large  acceleration  error.  With  the  accuracy  of 
the  angle  measurement  increased,  the  target  acceleration  estimate  after  the  maneuver 
onset  improves  faster  than  the  filter  that  uses  poor  angle  measurements.  This  is 
shown  in  Fig.  8.  However,  after  the  target  started  to  maneuver  in  the  middle  of 
the  engagement,  position  and  velocity  error  estimates  do  not  reduce  as  Vi  becomes 
smaller. 

Miss  distances  have  been  calculated  on  the  basis  of  50  runs  of  Monte  Carlo  sim¬ 
ulations  with  an  approximate  error  ±0.02  ft  due  to  subdiscretization  near  the  final 
time.  In  the  Table  1,  the  actual  states  are  fed  to  the  guidance  law  in  the  Case  I, 
and  the  estimated  states  and  maneuver  rate  estimate  are  fed  to  the  guidance  law  in 
the  Case  II  and  III.  The  estimates  are  obtained  from  angle-only  measurements  in  the 
Case  II,  and  from  both  angle  and  pseudo-measurement  in  the  Case  III.  Miss  distance 
performance  is  tested  as  more  noise  is  introduced  into  the  measurement  and  then 
into  the  dynamics.  For  the  particular  scenario  chosen  here,  miss  distance  has  been 
improved  by  using  the  angle  and  pseudo-measurement,  especially  as  the  process  noise 
power  spectral  density  0  in  the  state  dependent  noise  term  decreases. 

7.  Conclusions 

The  orthogonality  between  the  target  acceleration  and  velocity  vectors  is  a  typical 
characteristic  of  the  target  of  an  air-to-air  missile,  and  it  is  utilized  in  the  develop¬ 
ment  of  a  new  stochastic  target  acceleration  model  for  the  homing  missile  problem. 
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In  addition,  this  characteristic  is  also  implemented  in  the  form  of  an  augmented 
pseudo-measurement.  A  guidance  law  that  minimizes  a  quadratic  performance  index 
subject  to  the  stochastic  engagement  dynamics  is  determined  in  closed  form  where 
the  gains  are  an  explicit  function  of  the  estimated  target  maneuver  rate  and  time  to 
go.  Preliminary  results  for  the  two-dimensional  case  indicates  that  the  circular  target 
model  is  able  to  produce  a  reliable  estimate  in  the  homing  missile  engagement.  When 
it  is  augmented  by  the  fictitious  measurement,  the  modified  gain  extended  Kalman 
filter  using  the  proposed  target  model  results  in  the  remarkable  enhancement  of  tar¬ 
get  state  estimation  for  both  a  maneuvering  and  nonmaneuvering  situation.  This  is 
because  the  circular  target  model  inherently  includes  a  better  approximation  to  the 
simulation  dynamics  of  target  motion,  and  the  pseudo-measurement  imposes  satisfac¬ 
tion  of  the  orthogonality  characteristic  through  the  measurement  process.  However, 
the  improvement  in  miss  distance  for  a  particular  engagement  shows  significant  im¬ 
provement  in  going  from  angle-only  measurement  to  both  angle  measurement  and 
pseudo-measurement  when  the  process  uncertainty  decreases. 
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Appendix  A 

2.  Linear  Quadratic  guidance  law  for  deterministic  circular  target  model 

In  the  following,  the  optimal  deterministic  guidance  law  for  linear  quadratic  prob¬ 
lem  is  sought  for  the  current  circular  target  model  filter.  The  deterministic  optimal 
solution  can  be  obtained  by  solving  the  Riccati  equation  without  the  A  term  via 
transition  matrix  approach,  but  the  use  of  Euler-Lagrange  equation  seems  simpler  for 
this  case. 

The  problem  is  to  minimize  the  performance  index 

J  =  +  y)\  +  2 1  +  (A-1) 

subject  to  the  following  linear  dynamic  system 


X 

— 

u 

y 

= 

V 

ti 

= 

&Tx  ~  OM, 

v 

= 

ary  -  aMv 

dr. 

= 

0 

~~2°Tx  -uarv 

drv 

= 

0 

~2a7V  +  war. 

This  linear  system  of  dynamics  stems  from  taking  Ito  derivative  of  the  corre¬ 
sponding  nonlinear  stochastic  target  model  (6).  The  u  is  the  angular  rate  of  target 
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maneuver  which,  is  handled  as  a  known  constant  in  the  derivations.  In  the  acturl 
mechanization  of  guidance  command  the  value  of  Cj  constructed  from  the  estimated 
states  are  used. 

The  variational  Hamilt'.nan  and  the  augmented  end-point  function  are  given  by 

H  =  ca^+ca^  +  Aiu  +  Aju  +  A3(ar, -aw,) 

0  0 

+^4(07,  ~  “**»)  +  —  war,)  +  ^e(— -sr  +  uaTx) 

a  =  j(*?  +  v))  (A-4) 

where  A,-  , i  =  1,...,6  is  a  Lagrange  multiplier.  The  Euler-Lagrange  equations  for  A, 
are 

Aj  =  0,  A2  =  0,  —A3  —  Aj,  —  A4  =  A2  ’  (A-5) 

where  the  optimal  control  satisfies  the  optimality  condition 

A3  ^4 

aM,  =  — »  aM,  =  —  (A-6) 

Finally,  the  Euler-Lagrange  equations  with  the  natural  boundary  conditions  yield 

At  =  A2  =  yj,  A3  =  XfTgo ,  A4  =  y/Tgo  (A-7) 

which  gives  the  consol 

=  x{tj)Tgo/c,  aM,  =  y(tf)Tg0/c  (A-8) 

where  Tgo  is  the  time-to-go  of  missile  to  intercept  the  target  and  c  is  the  guidance  law 
design  parameter.  In  order  to  get  the  guidance  law  in  terms  of  the  current  states,  the 
underlining  dynamics  is  integrated  backward  from  tj  to  t.  Succesive  integrations  of 
state  differential  equations  yield 

aTx  =  vcosuTg0eTT'°aTl(tj) +u>sinuTg0e7T*°aTg{tj) 
ar„  =  -uJsinuTg0e?T'oaTt(ts)  +  ucosuTgoeTT’0  aT)l{t  j) 

u  =  hT°°X{tj) + u(*;) 
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+  e~--— [f  (1  -  cosuTgoeTT'°)  - u>sinuTgoeTT9°]aT,(tf) 
s-  +u>2  l 

+  e*“7— ~  cosu}Tgoe*'v'°)  -  ^-sinwTgoe%T'a)\aTy(tj) 

\  +  w2  ^ 

V  =  ^TaW/)  +  v(*/) 

+-ST-  rMl  -  cosuTgoe7T*°)  -f  ^r$iTUjjTg0e?T'o]aTx(tj) 

T  +  w  1  (A-9) 

4-sr^— r[?(l  -  cosuTgoe%T'°)  -  usinuTgoe*T'0))aTy{t}) 

X  +  u 3  1 

x  =  (!  -  ^)*(</)  -  Tg0u{tj) 

+  ^T+tJ^~~2££  “  (i+w2)(1  “  cos“T^T'°)  +  (<^  +  J) ^nu;rsoeTT^]aTl ',*/) 
+  el~-r[wTso  +  --(1  -  cosuTgoe?T'°)  +  smu>TgoefT'°)}aTy(tj) 

»  =  (l-|fM</)- !>(</) 

-  (f^(1  ”  c°sur>°c%T”)  -  (a;^{sinu,r-'-efT,'1°I'-(1,) 

+  91^  -  ^——^(1  -  cosuT10e‘T'")  +  -  ^  ■-sinuT,„erT”\<iT>{tt) 

X  +  w  l  — +W2  —  +  W2 

The  final  states  being  expressed  in  terms  of  the  current  states  via  6x6  matrix 
inversion, the  optimal  guidance  law  is  obtained  as  equation  (53)-  As  expected  from  dy¬ 
namic  coupling  in  the  target  acceleration  model,  guidance  commands  in  ~ach  channel 
axe  the  function  of  acceleration  components  in  both  x  and  y  directions. 


Appendix  B 


1.  State  and  Error  variance  associated  with  target  acceleration 


Since  the  initial  values  for  the  state  estimates  ascociated  with  target  acceleration 
are  set  to  zero,  the  state  and  error  variances  axe  computed  with  the  aid  of  expected 
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values  of  trigonometric  functions  such  as 

i?[cos20]  =  [  cos2  9p(9)d9  with  p(Q)  =  ,  ^  (B-l) 


By  using  standard  manipulation,  the  expected  value  of  the  cos2  5  is 

£[cos2  0]  =  ^[1  +  cos2&T26t],  (B-2) 

and  in  the  same  manner 

Elsin’Dl  =  i[l  +  a»2Je-"‘], 

£[cos  0  sin  0]  =  isin20e-2e‘. 

This  yeilds  the  initial  conditions  for  the  state  and  error  variances  associated  with 
target  acceleration  as  follows. 

P55(0)  =  X5S(0)  =  4maa(l  +  cos29)/2, 

Psa(0)  =  ^56(0)  =  4mmss*n  2^/2, 

/VO)  =  x57(0)  =  4mmaQ(l  +  cos29)/2, 

Pss(0)  =  XS8(0)  =  4m„Qsin  2^/2, 

P59(0)  =  X59(0)  =  4maiw2(l  +coS20)/2, 

•fsio(O)  =  X5io(0)  =  OjmaMU)2  sin29 /2, 

/VO)  =  ^6e(0)  =  Ormax(l  —  cos29)/2 , 

^7(0)  =  ^7(0)  =  ajmaJosin20 12, 

/VO)  =  *«(0)  =  -  cos20)/2, 

P.(0)  =  X69(0)  =  u>2sin29/2 , 

P«io(0)  =  X61o(0)  =  aO(l  "  cos20/2,  (B-4) 

P 77(0)  =  Xjj(0)  =  ajmMw2(l  +  cos29)/2, 

Pts(0)  =  X78(0)  =  a^maiu2sin29/2, 

/>79(0)  =  AVg(O)  =  dj.mMw3(l  +  cos29)/2, 

Prio(O)  -  A'71o(0)  =  a2Tmau3sin29/2, 

Ps8(0)  =  ^88(0)  =  4m„w2(l  -cos29)/2, 

Ps9(0)  =:  ^89(0)  =  aj-matuj3sin29 /2, 

P8io( 0)  =  X81o(0)  =  -  cos29)/2 , 

PVO)  =  ^99(0)  —  aTm«*u;',(l 

/^9io(0)  =  Agio(O)  =  (4maiu)*sin29 /2, 

/>ioio(0)  =  Aioio(O)  =  a7>mMcj'4(l  —  cos29)/2, 
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Table  1 


Case 

Range M  (ft) 

Miss  Distance(ft) 

e  =  o.oi 

V1  =  V0x  10-2 
V2  =  106 

I 

ebmm 

6000 

0.35 

0.32 

II 

6000 

1.30 

K2XJMI 

4000 

0.74 

III 

BKHHHi 

6000 

0.82 

mmwmm 

4000 

0.54 

0  =  0.001 

V1  =  Vo 

V2  =  106 

II 

6000 

2.65 

HHi 

4000 

1.71 

III 

KMfll 

6000 

2.13 

ksh 

4000 

1.36 

0  =  0.01 

Vi  =  Vo 

y2  =  io6 

II 

6000 

4.59 

4000 

1.97 

III 

HKHHHi 

6000 

4.29 

6000 

4000 

1.82 

0  =  0.1 

VJ  =  Vo 

V2  =  106 

II 

iHEMiHi 

6000 

4.98 

4000 

2.07 

III 

■Mi 

6000 

4.86 

6000 

4000 

2.01 

0  =  0.001 

Vi  =  Vo 

V2  =  108 

II 

6000 

6000 

2.65 

6000 

4000 

1.71 

III 

6000 

6000 

2.34 

6000 

4000 

1.49 
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YAXIS 


MISSILE-TARGET  TRAJECTORY 


Engagement  1  ( Ri  =  6000 ft,  Rm  =  6000/f) 


Engagement  2  (Ri  =  6000 ft,  Rm  —  4000/tf) 


Fig.  3  Typical  Missile-Target  Trajectories 

Ri  =  Initial  range,  Rm  =  Maneuver  onset  range 
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Fig.  4  Circular  target  model  filter  with  different  0's 
Ri  =  6000 ft,  Rm  =  6000ft,  cj  =  0.3 


:o 


ERROfl(f T/SEC)  CMON(rT) 


1  meaiwtimnl  j  ) 

V*  *  10*  |  I  . 

1/  -  .Ai  1*U  .  I  **  « 


— J - 1 - 1 - 1 _ lJ 

1-0  1.3  2.0  2.3  3.0 

TIlC  (  SCC) 


VELOCITY  ERROR 


TARGET  OMEGA  ESTIMATE 


Circular  target  model  filter  with  pseudo-measurement 
Ri  =  6000 ft,  Rm  =  6000/f,  u  =  0.3 
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Circular  target  model  filter  with  pseudo-measurement 
(A  case  with  different  V[s ) 

Ri  =  6000 ft,  Rm  =  4000/ 1,  u  =  0.3 


